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ABSTRACT 

We present results from a 49.3 ks Chandra X-ray observation of the strongly interacting dumbbell galaxies 
NGC 4782(3C278) and NGC 4783 that constrain the kinematics of the interaction and models for bending the 
radio jets associated with NGC 4782. The galaxies are embedded in an approximately spherical distribution of 
group gas, centered on NGC 4782, that, away from the interaction region, is nearly isothermal with mean tem- 
perature kT ~ 1.4±0.4keV. The X-ray morphology suggests that NGC 4783 is infalling into a single, massive 
galaxy group (LGG 316) with NGC 4782 nearly at rest at the center of the group potential. NGC 4783 shows 
a sharp X-ray surface brightness edge (cold front) to the east and a ~ 15 kpc ram-pressure-stripped tail to the 
west. Analysis of this cold front indicates NGC 4783 is moving east with a total velocity 870+4,™ km s" 1 (Mach 
1 -4^o 7) at an inclination angle 46°(> 33°) towards us with respect to the plane of the sky. A ~ 45 Myr old 
X-ray cavity, with enthalpy of 4.4 x 10 57 ergs, coincides with the eastern radio lobe of 3C278. X-ray knots are 
found on both the radio jet and counter-jet, coincident with peaks in the radio emission. If we assume a light, 
mildly relativistic jet in 3C278, then ram pressure velocities of 100- 200km s -1 impacting the eastern jet and 
~ 170 km s -1 acting on the western radio lobe are sufficient to produce their observed bending. These velocities 
may be caused by bulk motions of the surrounding gas induced by the high velocity interaction between the 
galaxies, by the motion of the host galaxy NGC 4782 relative to the IGM, or a combination of these processes. 
Subject headings: galaxies: groups: general - galaxies:individual (NGC 4782, NGC 4783) - galaxies: inter- 
galactic medium - X-rays: galaxies 



1. INTRODUCTION 

Deep, high angular resolution X-ray observations with 
Chandra and XMM-Newton have shown that the signatures 
of gas-dynamical interactions, i.e. shocks, cold fronts, wakes 
or tails, provide powerful new tools to constrain the three di- 
mensional motion of the galaxies with respect to the surround- 
ing IGM (Merrifield 1998; Vikhlinin et al. 2001), as well as 
revealing the thermodynamic state of gas in and near the in- 
teracting system. In particular, a leading cold front - trail- 
ing wake morphology of a galaxy moving through the sur- 
rounding IGM fixes the direction of motion of the galaxy in 
the plane of the sky. Analysis of cold fronts, first applied 
to merging subcluster cores by Vikhlinin et al. (2001), de- 
termine the three-dimensional velocity of galaxies infalling 
into nearby groups and clusters (e.g. Machacek et al. 2005a, 
Scharf et al. 2005 for NGC 1404 in the Fornax Cluster; 
Machacek et al. 2006 for NGC 4552 in the Virgo Cluster; 
Sun & Vikhlinin 2005 for UGC 6697 in the Abell 1367 clus- 
ter; Sun et al. 2005 for NGC 1265 in the Perseus Cluster). 
The morphology and thermodynamic properties of X-ray gas 
in the galaxy wakes and tails have been used to identify high 
velocity interactions and constrain the three dimensional mo- 
tion of the galaxy (Machacek et al. 2005b for NGC 6872 in 
the Pavo Group), as have the properties of shocks (for a re- 
cent example, see Rasmussen et al. 2006 for NGC 2276 in 
the NGC 2300 group). In interacting galaxies where gravi- 
tational tidal effects are also important, these gas-dynamical 
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methods provide important additional constraints on the dy- 
namics of the system that complement those obtained from 
modeling the tidal stellar distortions. 

Dumbbell galaxy systems (Wirth et al. 1982) can rep- 
resent high velocity, deeply inter-penetrating encounters of 
two roughly equal mass elliptical galaxies, with characteris- 
tic projected separations of ~ lOkpc (Prugniel & Davoust, 
1990) and mean relative radial velocity differences of ~ 
650km s" 1 (Valentijn & Casertano, 1990). In addition, dumb- 
bells frequently host radio sources in one or both of the galax- 
ies (e.g. Harris & Roberts 1960; Wirth et al. 1982; Parma 
et al. 1991; Gregorini et al. 1994). The majority of radio- 
active dumbbells are found near the cores of rich, often dy- 
namically young, clusters (see, e.g. Beers et al. 1992 for 
Abell 400; Reid et al. 1998 for Abell 3528). A few also 
have been found in galaxy groups. Since groups, like galaxy 
clusters, harbor dense X-ray gas (e.g. Schwartz, Schwarz & 
Tucker 1980; Mulchaey & Zabludoff 1998; Osmond & Pon- 
man 2004), hydrodynamical interactions between the galax- 
ies, radio jets, and group IGM in these dumbbells can be sig- 
nificant. 

In this work we present the results from a Chandra obser- 
vation of the radio-active interacting dumbbell system con- 
taining elliptical galaxies NGC 4782 and NGC 4783, lo- 
cated near the center of the galaxy group LGG 316 (Garcia 
1993), with line-of-sight velocities v r = 4618 ± 31 kms" 1 and 
v,. = 3995 ±22 kms" 1 , respectively (Martimbeau & Huchra, 
2006). Optical photometry shows that the central isophotes in 
each of the dumbbell galaxies are off-centered and the stellar 
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velocity profiles distorted, signaling the effects of intense tidal 
interactions between the pair (Borne et al. 1988, Madejsky 
et al. 1991, Madejsky 1992, Combes et al. 1995). Using dif- 
ferent numerical methods, Borne et al. (1988) and Madejsky 
& Bien (1993) obtained orbital solutions for NGC 4782 and 
NGC 4783 by fitting the observed isophotal shapes and stellar 
velocity distributions to simulations of tidal forces acting on 
the galaxies during the encounter. However, these solutions 
were not unique and, although these authors agreed that the 
data were consistent with a high relative velocity interaction 
viewed shortly after pericenter passage of the galaxies, they 
disagreed on the details of the orbital solution for the pair. 

NGC 4782 harbors the FRI radio source 3C278 (Harris & 
Roberts 1960) with two moderately bent radio jets and lobes 
(Baum et al. 1988) and a 5 GHz flux density for its point- 
like nucleus of 82 mJy. In contrast, VLA measurements give 
an upper limit on the 5 GHz flux density from the nucleus of 
NGC 4783, its dumbbell companion, of < 0.5 mJy (Hardcas- 
tle et al. in preparation). Early simulations of the evolution of 
NGC 4782's radio jets assumed the orbital solution of Borne 
et al. (1988), based on modeling the observed stellar tidal dis- 
tortions of the dumbbell assuming two approximately equal 
mass field galaxies initially infalling from infinity. Hot ISM 
gas distributions, associated with each of the galaxies, were 
introduced as free parameters in these simulations to provide 
ram pressure to account for the observed bending of the radio 
jets (Borne & Colina 1993). X-ray observations allow a di- 
rect measurement of the density distribution of hot gas in and 
surrounding each galaxy, reducing the number of free parame- 
ters needed by simulations. X-ray observations may also con- 
strain the motions of the interacting galaxies, allowing a more 
realistic choice of initial conditions for simulating the interac- 
tion. The presence of hot X-ray gas in NGC 4782/NGC 4783 
was first identified using the Einstein Observatory (Roberts 
et al. 1991; Fabbiano et al. 1992). Later ROSAT HRI ob- 
servations revealed a highly asymmetric hot gas distribution, 
with extended emission to the west of NGC 4783 and to the 
east of NGC 4782, and a bright bridge of X-ray emission con- 
necting the two galaxy centers (Colina & Borne 1995). In this 
work we use Chandra observations of NGC 4782/NGC 4783 
to measure the X-ray gas properties both in the interacting 
galaxies and the surrounding group, in order to constrain the 
dynamics of the system and the evolution of the radio features. 

This paper is organized as follows. In ^2] we discuss the 
observation and our data analysis procedures. In §|3] we ex- 
amine the effects of the high velocity encounter of NGC 4783 
with NGC 4782. We first measure the thermodynamic prop- 
erties of hot gas in the inner regions of the system (j j3. It and 
then characterize the mass, density and temperature structure 
of the larger scale group (j j3.3l l. In §|4] we discuss the X-ray 
features associated with the radio jets. In §|5] we use these 
data to constrain the kinematics of the galaxy interactions 
with each other and with the IGM, and present a physical pic- 
ture for NGC 4782/4783 dumbbell interaction suggested by 
our analysis in §|6] We summarize our primary findings in 
$7] Throughout this paper we adopt a flat, dark energy dom- 
inated cosmology (fi m = 0.3, A = 0.7) with Hubble parame- 
ter H = 70 km s" 1 Mpc" 1 . Taking the redshift of NGC 4782 
(z = 0.0154, Martimbeau & Huchra 2006) as representative of 
the group, the luminosity distance to the system is 66.7 Mpc 
and 1" corresponds to a distance scale of 0.314kpc. Unless 
otherwise indicated, all errors correspond to 90% confidence 
levels and coordinates are J2000.0. 



2. OBSERVATIONS AND DATA REDUCTION 

2.1. Chandra Observations 

Chandra observed the dumbbell galaxies NGC 4782 and 
NGC 4783 (obsid 3220) on 2002 June 16 for 49.3 ks using 
the Advanced CCD Imaging Spectrometer (ACIS, Garmire 
et al. 1992) in VFaint mode with the back-illuminated chip 
S3 at the aimpoint. The data were reprocessed using stan- 
dard X-ray analysis packages (CIAO 3.2, FTOOLS). Events 
were filtered to remove bad grades (1,5, 7), and hot pixels 
and columns. Additionally, we used VFaint mode filtering to 
obtain a factor ~ 2-3 better rejection of particle backgrounds 
at energies below 1 keV. The data were reprocessed using up- 
dated gain tables and corrected for the slow secular drift of the 
average pulse-height amplitude (PHA) values for photons of 
fixed energy (tgain 1 ) and for the buildup of contaminants on 
the optical filter (Plucinsky et al. 2003). Periods of anoma- 
lously high background (flares), as well as those with anoma- 
lously low count rates, were removed using the script lc_clean 
on the S3 data in the 2.5-7keV energy band after masking 
out bright sources. This resulted in a useful exposure time of 
46,667 s. 

We created backgrounds for imaging and spectral analyses 
of the group IGM from the 450 ks period D source free dataset 
aciss_D_7_bg_evt_010205.fits appropriate for the date of ob- 
servation and instrument configuration 2 . We checked the nor- 
malization of the source free background by comparing count 
rates in the source and background files in the 9.0- 1 1.5keV 
energy bands, where particle background dominates. We 
found that these count rates agreed to better than 1%, such 
that no additional background renormalization was required. 

We used a multiscale wavelet decomposition algorithm 
with a 5(7 detection threshold to identify X-ray point sources 
in the 8'. 4 x 8'.4 field of view of the S3 chip in four en- 
ergy bands: 0.3-1 keV (soft), 1 - 2keV (medium), 2-8keV 
(hard), and 0.3-8 keV (broad). We detected 59 X-ray sources 
on the S3 chip, in addition to the nucleus of NGC 4782, with 
29 of them lying within a 2' circle centered on NGC 4782. 
Three of these sources (two in the western jet and one in 
the eastern jet) are X-ray knots coincident with peaks in the 
radio emission and are discussed in 94.21 The bright point 
source at (12 A 54 m 35.4U ,-12°34'5.38"), located just north of 
NGC 4782's nucleus and containing 51 ±7 counts (la uncer- 
tainties) in the 0.3 - 8 keV energy band, is likely a ULX in 
the core of NGC 4782. For power law spectral photon indices 
r ~ 1 .6 - 2 and Galactic absorption, we find a 0.3 - 8 keV lu- 
minosity of (3 ± 1) x 10 39 ergs for this ULX. All point sources, 
except for the nucleus of NGC 4782, were removed from the 
surface brightness and spectral analyses. 

2.2. VLA Data 

To complement the Chandra data, we obtained previously 
unpublished deep observations from the archive of the NRAO 
Very Large Array (VLA), made at effective frequencies of 1 .5 
and 4.9 GHz. Details of these observations are given in Table 

m 

The VLA data were reduced in the standard manner within 
AIPS. Data from each individual observation were calibrated 
and self-calibrated; images made from long-baseline data 
were then used to cross-calibrate the short-baseline data and 
the datasets were concatenated to make a single dataset from 

1 see Vikhlinin et al. in http://cxc.harvard.edu/contrib/alexey/tgain/tgain.html 

2 see http://cxc.hai-vard.edu/contrib/maxim/acisbg 
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FIG. 1. — (top) 0.3 — 2 keV Chandra image of the dumbbell galaxies 
NGC 4782 and NGC 4783. The image has been background subtracted, 
exposure corrected and smoothed with a 1" Gaussian kernel. 1 pixel = 
0".49 X 0".49. (bottom) 03 -2keV Chandra image (as in the left panel) 
with optical contours superposed from a WFPC2 (FR680N) Hubble Space 
Telescope image. Contour levels correspond to flux densities of (3, 4.5, 6, 

o 

12, 24, 48) X 10" 21 ergs s _1 cirT 2 A~ l pixel -1 , respectively. 

which maps at various resolutions could be derived. In this 
paper we use the 1.5-GHz images to show the large-scale 
structure of the source, while the higher resolution available 
at 4.9 GHz is used to show details of the jet. 

The VLA data, together with our new Giant Meter-wave 
Radio Telescope (GMRT) observations of 3C 278, will be dis- 
cussed in more detail in a future paper. 

3. THE HIGH VELOCITY ENCOUNTER OF NGC 4782 
AND NGC 4783 

In Figure Q] we show the 0.3-2keV Chandra X-ray image 
of the interacting galaxies NGC 4782 and NGC 4783. The X- 
ray image has been background subtracted, corrected for tele- 
scope vignetting and detector response, and then smoothed 
with a 1" Gaussian kernel. In the bottom panel of Figure Q] 
we show isophotes from a WFPC2 (FR680N) Hubble Space 
Telescope image of NGC 4782 and NGC 4783 superposed on 
the 0.3-2keV Chandra image. 

We see that the peaks of the X-ray emission coincide with 
the off-centered peaks in the optical light at the locations of 
the galaxies' nuclei. Off-centered stellar isophotes are the 
defining signature of strong tidal interactions between ellip- 
tical galaxies, when the distance of approach is smaller than 
the characteristic size of the system. The off-centering results 
from the tidal deceleration of the outer stellar envelopes of 
each galaxy towards its partner, while stars near the galaxy 
centers have had time, due to their short (~ 10 7 yr) dynamical 



time scales, to average out the tidal effects and return to equi- 
librium (Madejsky et al. 1991, Combes et al. 1995). More 
asymmetric extended features, such as tidal tails, are not yet 
seen, since the stars which became unbound during the peri- 
center passage of the galaxies, have not had time to escape the 
parent galaxies (Combes et al. 1995). 

A key feature that distinguishes gravitational interac- 
tions (tides) from gas-dynamical processes (ram-pressure or 
turbulent-viscous stripping) is that gravitational interactions 
act on both stars and gas, while gas-dynamical processes 
are ineffective at distorting the stellar distributions or veloc- 
ities. From the bottom panel of Figure Q] we see an en- 
hanced 'bridge' of X-ray emission, along the line (oriented 
19° east of north) connecting the galaxy centers. The 'bridge' 
is closely correlated with the stellar distribution, as shown 
by the HST isophotes (overlaid contour lines), and coincides 
with the direction of the north-south tidal distortions in the 
stellar isophotes. This suggests that the formation of the X- 
ray 'bridge' is due in part to the same tidal interactions that 
disturbed the stars. Emission in the 'bridge' region is also 
likely partly due to projection effects from the superposition 
of the two galaxy ISMs caused by the closeness of the en- 
counter. 

The X-ray surface brightness distribution has extended fea- 
tures that do not follow the stellar distribution, suggesting that 
hydrodynamical processes are important. In the top panel of 
FigureQ] NGC 4783 shows a sharp surface brightness discon- 
tinuity (edge) on the eastern side of the galaxy, a ~ 4- 5 kpc 
long streamer or 'horn' of gas swept back from the edge to the 
northwest, and a tail of enhanced X-ray emission, extending 
~ 15 kpc to the west, opposite NGC 4783's 'leading' edge. 
These features are similar to those predicted in numerical sim- 
ulations of gas-dynamical stripping of elliptical galaxies due 
to their motion through the ambient IGM (Stevens etal. 1999; 
Toniazzo & Schindler 2001; Acreman et al. 2003; Heinz 
et al. 2003), and to features seen in other elliptical galaxies 
undergoing ram-pressure and/or turbulent-viscous stripping. 
The leading edge - trailing tail morphology determines that 
NGC 4783's motion in the plane of the sky is towards the east. 
This is in qualitative agreement with the relative galaxy mo- 
tions from the simulations of Borne et al. (1988) and Made- 
jsky et al. (1993), where NGC 4783 passes to the west of 
NGC 4782 at pericenter. The filamentary 'horn' extending 
from the sharp edge of NGC 4783 is likely gas being stripped 
at the galaxy-IGM interface. The twisted morphology of the 
bright inner regions of NGC 4783's tail is reminiscent of von 
Karman vortices seen in hydrodynamic fluid flows and pro- 
duced in recent simulations of ram pressure stripping of disk 
galaxies (Roediger et al. 2006). 

An asymmetric 'fan' of X-ray emission extends to the east- 
northeast of NGC 4782, overlapping in projection the bend in 
the eastern radio jet (see the bottom panel of Fig. |4j, while 
less X-ray emission is observed to the southwest. However, 
identifying the dominant processes responsible for the ori- 
gin of the fan is complex. Madejsky (1992) found a sharp 
increase in the stellar velocity dispersion ~ 15-20" from 
NGC 4782, in the base of the fan, indicating the strong in- 
fluence of tidal forces on matter (gas and stars) in the region. 
The compression wave produced by the high speed passage of 
the companion galaxy NGC4783 from the west to the north- 
east of NGC 4782 would displace the ISM outside the core of 
NGC 4782 to the east relative to NGC 4782's stellar distribu- 
tion, as is seen. Similar bulk motions, on larger scales, have 
been shown to occur in ICM gas in merging clusters (Asca- 



4 



n 




FIG. 2. — (top) 0.3 — 2keV background subtracted, exposure corrected 
and Gaussian smoothed Chandra image of the dumbbell galaxies NGC 4782 
and NGC 4783 with spectral regions for N4783's head (H) and tail (T), the 
interconnecting bridge (B), eastern fan (F), and NGC 4782's nucleus (N) 
and surrounding annulus (A) superposed. Local background regions (dashed 
lines) are also shown for region H (circles), A (annular sector), T(white box), 
B(right yellow box), and F(left yellow box), {bottom) Gas temperatures in the 
elliptical galaxies NGC 4782 and NGC 4783 for the spectral regions shown 
to the left and spectral models listed in TablefJ] 

sibar & Markevitch 2006). Any small velocity of NGC 4782 
relative to the group potential and ambient IGM, induced by 
the tidal interaction of the galaxy pair, also would contribute 
to the hydrodynamical displacement of the gas. Finally, as we 
show in 34.21 some of the X-ray emission might be associated 
with the eastern radio jet. 

To probe the evolution of NGC 4782, NGC 4783 and the 
group IGM, we next characterize the thermodynamic prop- 
erties of the hot gas in the galaxies and tail, bri dge a nd fan 
( j33.lt . as well as in the surrounding group IGM ( 33.31 >. 

3.1. Galaxy Gas Temperatures and Densities 

The spectral extraction regions used to determine the prop- 
erties of hot gas in the elliptical galaxies NGC 4782 and 
NGC 4783 and in the nearby extended X-ray features are 
listed in Table [2] and are shown in the top panel of Figure 
[2] We use local background regions with the same geome- 
tries and at the same mean radial distances as the regions of 
interest to subtract the contribution of the surrounding group 
IGM, as well as the particle, soft Galactic, and cosmic X-ray 



backgrounds, from the galaxy spectra for regions H, T, B, and 
F. For the annular region A in NGC 4782, we use an annular 
sector abutting region A to the south. These local background 
regions are also listed in Table [2] and delineated with dashed 
lines in the bottom panel of Figure [2] For the background of 
the nuclear region N of NGC 4782, we used a small annu- 
lus (region A in Fig. |2). Counts were grouped for spectral 
fitting with a minimum 20 counts per bin. All spectra were 
then fit using XSPEC 1 1.3 with instrumental response matri- 
ces produced using standard CIAO tools. We restricted our 
spectral fits to the 0.3-8 keV energy range for regions H, T, 
and F, 0.3-2keV for regions B and A, and 0.3-5keV for 
the NGC 4782 nucleus (region N). We discuss our spectral 
analysis for all regions outside the nucleus below, and for the 
nuclear region in 33.21 

We used single temperature APEC thermal plasma models 
(Smith et al. 2001) with solar abundance tables from Anders 
& Grevasse (1989) to model the hot ISM. Each model com- 
ponent was corrected for absorption using Wisconsin photo- 
electric cross-sections (Morrison & McCammon 1983). We 
found no evidence for increased absorption over the Galac- 
tic value, and so fixed the hydrogen absorbing column at 
3.58 x 10 20 cm" 2 (Dickey & Lockman 1990) for all model 
components. Since abundances in nearby elliptical galaxies 
are known to vary, we explore models with fixed metal abun- 
dance between 0.3- 1.2Z Q , where the lower bound is taken 
from the 90% CL lo wer limit on the metallicity of the sur- 
rounding IGM (see 33.31 and Table |9]l and the upper bound 
corresponds to the highest metallicities observed in compara- 
ble, moderate-X-ray-luminous elliptical galaxies in the sam- 
ple of Humphrey & Buote (2006). 

Our regions contain considerable stellar mass, such that the 
contribution from LMXBs may be significant. We correct 
for the presence of LMXBs following Gilfanov (2004), who 
showed that the X-ray luminosity function for LMXBs in el- 
liptical galaxies has a universal shape with its normalization 
proportional to the galaxy's stellar mass. We determine the 
contribution of LMXBs to the 0.5 - 8 keV X-ray luminosity 
from the stellar mass in each spectral region, determined from 
K-band photometry. The contribution of LMXBs to each X- 
ray spectrum is modeled by adding a 5 keV bremsstrahlung 
component to the spectral models (Kraft et al. 2001) with its 
normalization fixed by the expected LMXB unabsorbed X- 
ray flux. We checked that our results were unchanged if we 
instead chose to model the LMXB contribution to the spectra 
with a (T = 1 .6 ) power law. 

The results of our spectral fits and the intrinsic X-ray lumi- 
nosities for each region are given in Tables [3] and |4] respec- 
tively. Only in two regions (H and A) does the 0.5-2keV 
emission from LMXB's contribute more than 20% of the to- 
tal emission. We find that, once the contribution of LMXBs 
has been taken into account, single temperature APEC models 
for hot gas provide excellent fits to the remaining X-ray emis- 
sion in all regions outside the nuclear region (N). The best fit 
temperatures for these regions do not change when the abun- 
dances are varied between 0.3 - 1 .2Z Q . Similarly, the derived 
0.5-2 keV and total 0.5- lOkeV X-ray luminosities vary by 
<4%, andthe2-10keV luminosity by < 17%, over the abun- 
dance range of interest. Thus, for simplicity, we present only 
the 0.5Z Q metallicity model results in Tables [3] and [4] We 
plot the best fit gas temperatures as a function of region la- 
bel in the bottom panel of Figure |2] The gas temperature in 
the central regions of the galaxies (0.48+oJj keV for region 
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H in NGC 4783 and 0.43^ keV for the annular region A 
in NGC 4782) are typical of those found in the cores of el- 
liptical galaxies. Gas temperatures in the extended features, 
0.68^$ keV for NGC 4783 ' s tail re g ion T ' °- 68 -o:o8 keV for 
the bridge B, and 0.56^JJ [4 for the fan F are typical of galaxy 
gas that has been stripped by hydrodynamical or tidal pro- 
cesses and significantly cooler than the temperature measured 
for the group IGM (see $333\ . 

We use the best fit spectral models to estimate the mean 
density and gas mass in each region (see, e.g. Eq. 3 in 
Machacek et al. 2006), and to calculate the thermal gas pres- 
sures. For the tail region (T), bridge (B) and fan (F), we as- 
sume gas fills a cylinder of rotation about the major axis with 
(projected) length l p and diameter given in Table [2] We as- 
sume spherical volumes for the central regions of the galaxies. 
Our results for the electron densities, thermal pressures, cool- 
ing times, and total gas masses, assuming uniform filling and 
zero inclination angle with respect to the plane of the sky, are 
listed in Table [5] where the uncertainties reflect both the 90% 
CL uncertainties in the temperature from the spectral fits and 
the uncertainties in the assumed metal abundance. The upper 
and lower uncertainties on the density, pressure and gas mass 
are derived from the low (0.3Z Q ) and high (1 .2Z Q ) abundance 
models, respectively. The mean electron density in each re- 
gion also depends upon the relative geometry of the collision 
through the inclination angle Q of each region with respect to 
the plane of the sky. The physical length (and thus volume) 
of each cylindrical region (T, B, and F) scales as cos _1 (C)> 
such that the gas density and pressure decrease with increas- 
ing angle of inclination, scaling as cos 1/2 (Q, while the total 
gas mass increases, scaling as cos -1 /2 (Q. 

For zero inclination angle and uniform filling, we find mean 
electron densities n e ~ 0.01 -0.03 cm" 3 , gas pressures p ~ 
(1-5) x 10 _11 ergs cm" 3 , and cooling times ~ 0.2-1. 2 Gyr 
(long compared to the timescale of the collision), for all re- 
gions outside the nuclear region N. Hot gas masses in the 
central regions of the galaxies are - 2.9+% x 1O 7 M in the 
central 2kpc of NGC 4783 and - 4.4^;| x 10 7 M Q in the 
central 3.4 kpc of NGC 4782, comparable to gas masses 
(3.6 x 10 7 M Q for NGC 4783 and 8 x 10 7 M Q for NGC 4782) 
found by direct integration of the /3-model density distribu- 
tion fits of Jones et al. (2007). If we compare these cen- 
tral gas masses to those within NGC 4472, a noninteract- 
ing elliptical galaxy with similar velocity dispersion, we find 
that the central gas masses in NGC 4472 (i.e. 6.8 x 10 7 M Q 
and 1.43 x 10 8 M© within r = 2 and 3.4 kpc, respectively, of 
NGC 4472's nucleus; Jones et al. 2007) are a factor ~ 2- 3 
larger than the gas masses found in the central regions of the 
interacting galaxies NGC 4782 and NGC 4783, suggesting 
that gas may have been depleted in the central regions of the 
dumbbell galaxies as a result of the interaction. 

For nonuniform filling (77,- < 1, for each region z) the gas 
density and pressure are increased by a factor (fy) -1 / 2 and total 
mass reduced by (77,) '/ 2 . We caution the reader that in regions, 
such as region F, where the emission is clumpy and the filling 
factor likely < 1, the values in Table [5]represent lower bounds 
on the density and pressure and an upper bound on the derived 
X-ray gas mass in the region. 

3.2. The AGN at the Heart of NGC 4782 

Our spectral fits to the nuclear region N of NGC 4782, that 
hosts the radio source 3C278, are listed in Table|6] Neither a 
single temperature APEC model nor a single power law model 



can describe the spectrum. We obtained a good representation 
of the spectrum using an absorbed power law plus thermal 
APEC model. Due to our low statistics (193 ± 14.5 source 
counts), we fix the hydrogen absorption column at Galac- 
tic (3.58 x 10 20 cm" 2 ) and explore a range of photon indices 
(r = 1.2-2.2) and gas metallicities (0.3- 1.5Z Q ), typical for 
LLAGN (see, e.g. Evans et al. 2006 and references therein), 
and gas in the centers of elliptical galaxies (Humphrey & 
Buote 2006). Since the temperature of the gas is unchanged 
(at 90% CL) when the power law photon index (T) and gas 
metal abundance are varied over these ranges, we list only a 
single case (r = 1.7 and metal abundance 0.5 Zs) in Table [6] 
We find a temperature kT = 0.4+q 2 keV, signaling the presence 
of cool X-ray emitting gas within 0.7 kpc of NGC 4782's nu- 
cleus (See Fig. |2). The best fit temperature for the thermal 
component is also unchanged when the bremsstrahlung com- 
ponent, accounting for LMXB's, is added to the model. 

In the nuclear region of NGC 4782, the total 0.5 - lOkeV 
luminosity of (1.5 ±0.3) x 10 40 ergss _I is dominated by the 
power law component with (9 ±2) x 10 39 ergss _1 , 4.2^'| x 
10 39 ergs s" 1 and 2.1 x 10 39 ergs s _I attributed to the power 
law, thermal (gas), and bremsstrahlung (LMXB) components, 
respectively. Since power law X-ray emission may be pro- 
duced by subparsec scale jets in the core of radio sources, 
as well as by accretion onto the central black hole (Evans 
et al. 2006), the observed power law X-ray luminosity ( 
(9 ± 2) x 10 39 ergs s" 1 ) is an upper bound on the X-ray lu- 
minosity of the central engine. 

We estimate the mass (Mbh) of the central black hole in 
NGC 4782 to be 1O 9 M using the central stellar velocity dis- 
persion of 3 10km s" 1 (Madejsky et al. 1991; Madejsky 1992) 
in the Mbh - &c relation (Ferrarese & Merritt 2000; Gebhardt 
et al. 2000). The Eddington luminosity for such a massive 
black hole is LEdd ~ 1.3 x 10 47 ergs s" 1 , such that Lx/^Edd < 
10~ 7 . Even if the bolometric luminosity is 3-10 times higher 
than the observed 0.5- lOkeV X-ray luminosity, as is the case 
for accretion-dominated AGN (Elvis et al. 1994), the central 
AGN in NGC 4782 is radiating at very low luminosity for its 
mass. 

Variation of the metal abundance does, however, affect 
the APEC component normalization, which decreases as the 
model metallicity is increased, and thus the derived gas den- 
sity, thermal pressure and gas cooling times. These quantities 
remain insensitive to the variation of the power law photon in- 
dex. For all quantities below, the quoted uncertainites reflect 
both the uncertainties due to the variation of gas metal abun- 
dance and LLAGN photon index in the spectral model and the 
90% CL uncertainties in the gas temperature in the model fit. 
We find a mean electron density n e ~ 0.09 ± 0.04 cm" 3 and 
thermal pressure p ~ (1.2^ 7) X 10~ 10 ergs cm" 3 . The cool- 
ing time (~ 80^ Myr) in region N is short, suggesting that 
gas accretion onto the supermassive black hole can trigger 
episodes of AGN jet activity. 

Given the density and temperature of the gas in the nuclear 
region, we can estimate the power available from Bondi ac- 
cretion (Mbhc 2 ) onto the central black hole, where c is the 
speed of light. The Bondi accretion rate is 

Mbh = 7 x 10 23 

x (M BH /lO 9 M ) 2 (n e ,/O.17cm- 3 )(A:r/O.8keVr 3/2 gs- 1 (1) 

(di Matteo et al. 2003). Using the density (0.09 ±0.04 cm" 3 ) 
and temperature (0.4±Q 2 keV) from region N, we find a Bondi 
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accretion rate of M B h > OMO^^MQ-yr' 1 . This is a lower 
bound on the Bondi accretion rate, because the Bondi accre- 
tion radius (r a = GMbh/c 2 ) is much smaller than the 0.7 kpc 
radius of region N, such that the mean density from our spec- 
tral fit to region N likely underestimates the gas density at the 
accretion radius. Assuming the lower bound on the Bondi ac- 
cretion rate (0.005M Q yr ) from the lower limit on the gas 
density (metal abundance 1 .5Z ) and upper 90% confidence 
limit on the temperature, the minimum total power available 
from Bondi accretion (Mbhc 2 ) is then >3x 10 44 ergs s" 1 . As- 
suming the standard radiative efficiency of 0.1, we would ex- 
pect an accretion luminosity onto the central black hole in 
NGC 4782 greater than 3 x 10 43 ergs s" 1 . The observed X-ray 
luminosity is more than 3000 times smaller. The inferred ra- 
diative efficiency for the central AGN, L x /Mbhc 2 < 3 x 10" 5 , 
is then an upper bound to the true radiative efficiency, both be- 
cause the Bondi accretion rate may be underestimated and be- 
cause the observed nuclear X-ray luminosity for FRI sources 
is likely contaminated by emission from a subparsec scale jet 
(Evans et al. 2006). 

3.3. Characterization of the Group Environment 

NGC 4782 and NGC 4783 are not an isolated galaxy pair, 
but are part of the LGG 316 galaxy group (de Souza & Quin- 
tana 1990; Quintana et al. 1996). The derived properties of 
LGG316 are summarized in Tableland discussed below. 

3.3.1. Group Velocity Dispersion and Dynamical Mass 

We adopt a radial velocity of 4618 ± 31 km s" 1 for 
NGC 4782 and 3995 ± 22 for the northern companion 
NGC 4783 (Martimbeau & Huchra 2006), in agreement with 
Borne et al. (1988), Madejsky et al. (1992), and Quin- 
tana et al. (1996). From 22 galaxies with radial veloci- 
ties between 3500 km s -1 and 5500 km s -1 in a 1° square re- 
gion centered on NGC 4782, we found a mean radial ve- 
locity of 4604 km s" 1 and line-of-sight velocity dispersion of 
447 ±70 km s" 1 . Then using the mass-velocity dispersion 
(M - a) relationship normalized by simulations (f a ~ 0.9, 
Bryan & Norman 1998), we estimate a dynamical mass within 
r5oo for the group of ~ 7 x 1O 13 M with an uncertainty of 
~ 30%. These results are in excellent agreement with previ- 
ous results (de Souza & Quintana 1990; Quintana et al. 1996) 
when the galaxies are considered members of a single galaxy 
group (LGG 316) with NGC 4782 at its center. However, 
Quintana et al. (1996) also suggest that their data are con- 
sistent with a double-peaked galaxy distribution, in which 
NGC 4782 and NGC 4783 are the dominant members of two 
unequal mass subgroups undergoing a high velocity collision 
and possible merger. 

3.3.2. Surface Brightness Profile, Temperature and Density 

The characterization of the hydrodynamical interactions, 
acting to transform the galaxies during their encounter, re- 
quires a detailed understanding of the thermodynamic proper- 
ties of the surrounding group IGM. Furthermore, X-ray mea- 
surements of the hot gas density and temperature distributions 
in the group IGM may reliably trace the underlying dark mat- 
ter potential and reflect the dynamical state of the system. 
Large scale asymmetries in the temperature and density pro- 
files in the hot gas in the cores of galaxy clusters, such as those 
from cold fronts and shocks, have been used to identify super- 
sonic mergers of subcluster components in those systems (e.g. 
see Markevitch et al. 2002, 2005). Such signatures would be 



expected in the hot IGM of groups, such as LGG 316, if the 
ongoing encounter is the supersonic collision of two galaxy 
subgroups. 

In the upper left panel of Figure [3] we show the 0.3-2 keV 
image of the full ACIS S3 CCD. The galaxies are embedded in 
a roughly circular region of dense gas in the group's core, ex- 
tending ~ 50 kpc from the center of the dumbbell system. We 
find no visual evidence of asymmetries in the surface bright- 
ness outside the central interaction region, i.e. at r > 17 kpc, 
that might be expected for a high velocity merger of two sub- 
groups, as suggested by Quintana et al. (1996). 

In the upper right panel of Figure [3] we show the surface 
brightness profile as a function of radial distance from the cen- 
ter of NGC 4782 taken in two angular sectors (chosen to ex- 
clude regions containing the eastern and western radio lobes), 
with angles, measured counterclockwise from west, ranging 
from 38° to 88° for the northern profile (open squares) and 
from 200° to 336° for the southern profile (filled circles). 
Since ROSAT All Sky Survey maps of the Galaxy show excess 
soft X-ray emission near NGC 4782/4783, we construct the 
radial surface brightness profile in the 1-2 keV energy band 
where this soft Galactic background is negligible. To probe 
the group emission at larger radii, we also show the surface 
brightness determined from a 443" x 443" rectangular region 
on ccd SI (filled triangle). Outside the cores of the galaxies 
and regions of enhanced emission due to stripped galaxy gas, 
i.e. radii r > 15 kpc), we find that the surface brightness pro- 
files both to the north and south of the interacting galaxy pair 
can be described by a single power law, consistent with gas 
residing in the gravitational potential of a single, large group. 
From the surface brightness we infer an electron density for 
the IGM of niGM ~ 2.5^ x 1 0" 3 cm" 3 (r/ 15.45 kpc)- 3/3 with 
j3 = 0.35, where the uncertainties reflect the uncertainties in 
the spectral models for the IGM gas, determined by varying 
the gas temperature between 1.0-1. 8 keV and metallicities 
between 0.3 - 0.5Zq (see the lower left panel of Fig. [3). 

We first determined the mean temperature and abundance of 
the group IGM as a function of radius by fitting the spectrum 
in five annular regions (P„ P, P„i, P„2, P W 3; see the upper left 
panel of Fig. [3]and Tabled, covering as much of CCD S3 as 
possible (r < 100 kpc) outside the interaction region contam- 
inated by galaxy gas, and, unless otherwise noted, our discus- 
sion of thermal properties of the group IGM are based on the 
spectral fits to these regions. We found excess emission in the 
0.5-0.6keV band in each spectrum, consistent with a back- 
ground due to local emission from Galactic oxygen above 
that subtracted by the blank sky backgrounds. We modeled 
this soft excess background contribution to each spectrum 
by adding an additional local (z = 0) low temperature APEC 
component, fixing all of the APEC parameters for the soft 
Galactic oxygen component at the best fit values from a two 
component APEC fit to the spectrum of annulus P„i. Region 
P„i was chosen because it has good statistics (2269 ±67 net 
counts) and, with an inner radius of 37 kpc, is well outside the 
dumbbell galaxies, avoiding contamination by galaxy gas. We 
found a best fit temperature of 0. 1 8±q keV and metallicity of 
1 .2(> 0.2) Zq for the soft Galactic oxygen background at z = 
0, and a temperature and metal abundance of 1 .87 ± 0.28 ke V 
and O.28!g-gZ (x 2 / dof = 121/133) , respectively, for the 
LGG 316 group IGM (z = 0.0154) in that region. We found 
a mean 0.3-2keV surface brightness from the soft Galactic 
component of 5.7 x 10~ I9 ergss _I cm -2 arcsec -2 . 

We checked that this model for the local soft X-ray compo- 
nent was robust in two ways. First, since the determination of 
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FIG. 3. — (upper left) 0.3-2keV background subtracted, exposure corrected Chandra image of the region (chip S3) surrounding the dumbbell galaxies 
NGC 4782 and NGC 4783 with spectral regions superposed. The image has been Gaussian smoothed with a <r = 4" kernel to highlight the X-ray emission 
from the group IGM. The spectral regions are 40° annular sections numbered counterclockwise starting from 345°. Labels 4 and 1 mark the locations of the 
eastern and western radio lobes, respectively, (upper right) 1 — 2 keV surface brightness profiles as a function of the radial distance r measured from the center 
of NGC 4782 in an angular wedge from 38° to 88° to the north (open squares) and from 200° to 336° to the south (filled circles), measured counterclockwise 
from west. The filled triangle denotes the mean surface brightness in a rectangular region on ccd SI spanning radii 225 < r < 370 kpc (horizontal line). The 
diagonal line shows the surface brightness profile expected for gas with temperature (abundance) of 1.4 keV (OAZq), respectively, and electron density of the 
form n e = 0.0025(r/15.45kpc)~' 05 . (lower left) IGM temperature (circles) and abundance (triangles), with 90% confidence level uncertainties, as a function of 
radial distance r from the center of NGC 4782. (lower right) Gas temperature measured in the 40° sections for region P, shown in the upper left panel (see also 
Table|2), as a function of region label. Solid and dot-dashed lines denote the temperature from the simultaneous fit to the western segment (sectors 1 - 3 & 9) and 
eastern segment (sectors 4—8), respectively, with the 90% confidence levels for the western (eastern) segment shown as thin dashed (dotted) lines. 
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the abundance for the local emission is uncertain, we varied 
the abundance value for the local emission in region P„i from 
0.2 -2Z and found that the temperature of the soft Galactic 
component was unchanged (at 90% confidence limits) and the 
0.3-2keV fluxes for the soft Galactic component varied by 
< 10%. Second, assuming that the mean flux per pixel in the 
soft Galactic component is constant across the detector, we 
used our model based on region P„i to predict the 0.3 - 2 keV 
soft Galactic flux expected from a 443" x 443" rectangular 
region on ccd SI, located ~ 15'. 6 from NGC 4782, where 
the ratio of soft Galactic emission to group IGM emission is 
^2.3. We then fit the spectrum of this rectangular region with 
a two component APEC model, fixing the temperature and 
abundance (0.18keV and 1.2Z Q , respectively) from our P„i 
model for the soft X-ray background, but allowing the model 
normalization to vary. Metal abundances of the second APEC 
(IGM) component were varied between 0.1 -0.3Z Q , with the 
temperature and normalization determined from the fit. In all 
cases, the measured 0.3-2keV soft Galactic foreground flux 
from this region on ccd SI differed by less than 6% from that 
predicted by our background model rescaled from region P„i . 
Similarly, the temperatures found for the second APEC com- 
ponent (~ 1 -2keV) were consistent with emission from the 
group IGM. 

Thus we fix the normalization of the soft Galactic back- 
ground component for all other spectral regions at the best 
fit value for region P„i scaled by the ratio of the region ar- 
eas. We then fit the spectra for the annular regions, allow- 
ing the temperature, abundance and normalization of the IGM 
APEC component to vary. Our results are listed in Table [9] 
In the lower left panel of Figure [3] we plot the mean tem- 
perature (circles, solid line) and mean abundance (triangles, 
dashed line) from these model fits as a function of radial dis- 
tance from NGC 4782. Within the uncertainties, the group gas 
is nearly isothermal at ~ 1.4±0.4keV. The measured metal 
abundances of ~ 0.3-0.5Z Q are constant at the 90% confi- 
dence level. 

To investigate any angular change in the temperature of the 
surrounding group IGM, we divide each annulus into sectors 
(see the upper left panel of Fig. |3j, based on nine 40° sec- 
tors of the complete reference annulus P. The annular sectors 
in region P were chosen such that the western radio lobe lies 
mainly in annular sector 1 and the eastern radio lobe in annu- 
lar sector 4. For annulus P, we fit an absorbed single tempera- 
ture APEC model to the IGM emission in the individual annu- 
lar sectors, fixing the abundance at the best fit value (0.46Z Q ) 
for the annulus as a whole. Our results are plotted in the lower 
right panel of Figure [3] The data suggest that the IGM gas 
temperature may be higher to the west than to the east in this 
annulus (17 ^ r ^ 37kpc) . To test this we divided the an- 
nulus P into two segments, the eastern segment (annular sec- 
tors 4-8) and the western segment (annular sectors 1-3 and 
9), and simultaneously fit the spectra of the annular sectors 
in each segment, using an absorbed APEC model as above. 
Our results, shown as the dot-dashed and solid lines for the 
eastern and western segments, respectively, in the lower right 
panel of Figure [3] show that, at the 90% confidence level, the 
mean temperature of the gas in the western segment is higher 
(kT = 2.5 1±038 keV) than the IGM gas temperature in the east- 
ern segment at the same radii (kT = 1.75^q 13 keV). From fit- 
ting the spectra of the angular sectors in the remaining outer 
annuli , we find no evidence (at 90% confidence) for angu- 
lar variation in the IGM gas temperature at these larger radii 



(r > 37kpc). Thus the overall gas morphology and tempera- 
ture of the IGM for radii r>37 kpc is consistent with a sin- 
gle group of galaxies, in which the group gas at large radii is 
expected to be close to isothermal and in hydrostatic equilib- 
rium; while the increased temperature to the west, close to the 
interaction region, may be evidence for the recent supersonic 
passage of NGC 4783 to the west of NGC 4782, in agreement 
with simulations (Borne et al. 1988; Madejsky & Bien 1993), 
thereby heating the IGM. 

3.3.3. Mass, Luminosity and Group Scaling Laws 

We use the mean temperature of the group gas in the outer 
annuli in the M-T relation from Finoguenov et al. (2001) to 
estimate the mass of the group, 

M 500 = (2.64^ x 10 I3 )(£L) L78±(U . (2) 

From Equation [2] we find a group mass of 4. 8^2 x 1O 13 M 
for a mean temperature kT = 1.4±0.4keV, in agreement 
with the group mass determined from the (447 km s" 1 ) ve- 
locity dispersion for the LGG 316 group taken as a whole. 
The total X-ray luminosity of the group is more difficult to 
determine, since we can only directly measure the surface 
brightness profile and verify our model for the gas density 
to radii of ~ 50 kpc to the south and ~ 100 kpc to the north 
of NGC 4782, i.e. 0.08r 50 o and 0.15r 50 o, respectively. We 
determine the bolometric X-ray luminosity within radius r 
by using the 1-2 keV flux, found from integrating the (3- 
model for the surface brightness profile in Fig. [5J to com- 
pute the unabsorbed 0.01-14keV flux, using APEC spec- 
tral models with gas temperatures 1.4±0.4keV and abun- 
dances 0.3 ± O.2Z from our spectral fits to the IGM gas. 
We find the X-ray luminosity of the group gas within a ra- 
dius of 50(100) kpc of NGC 4782 to be logL x = 41. 7 1^04 
(41.99^J ( 1 ) 4)ergss" 1 , where the uncertainties in the luminos- 
ity reflect the uncertainties in the spectral model temperatures 
and abundances. However, as the upper right panel of Fig- 
ure [3] shows, data suggest that the (3 = 0.35 model for the 
X-ray surface brightness (and gas density) remains valid to 
r > 300 kpc (> 0.45r5oo), implying a lower bound on the X- 
ray luminosity for the group of logLx > 42.42±g 'J ergs s" 1 . 
Although the extrapolation of the surface brightness beyond 
r ~ 300 kpc may introduce large uncertainties, we use this 
extrapolation to estimate the total X-ray luminosity expected 
within f-500 = 660 kpc for comparison with X-ray scaling re- 
lations. We find logLx = 42.73^;^| ergs s" 1 for the LGG316 
group. This result is in excellent agreement with that expected 
from the Lx-7x and Lx-(T V scaling relations for a single large 
X-ray group of temperature 1 .4 ± 0.4 keV and velocity disper- 
sion ~ 447 km s" 1 (see Osmond & Ponman 2004, Figs. 13 and 
15, respectively). 

4. X-RAY FEATURES ASSOCIATED WITH THE RADIO 
JETS AND LOBES 

In the top panel of Figure HI we superpose contours from a 
1 . 5 GHz VL A map on the . 3 - 2 ke V Chandra image to show 
the correspondence between the large scale radio structure of 
3C278 and X-ray emission from NGC 4782/4783 and the 
LGG316 group. The appearance of the eastern and western 
radio jets and lobes are very different. The eastern jet shows a 
prominent bend towards the northeast at ~ 5-6kpc (17-18") 
from the nucleus of NGC 4782. Shortly after the bend, the 
eastern jet inflates into an approximately spherical lobe of ra- 
dius ~ 10.7 kpc (34"). The western jet extends straight for 
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FIG. 4. — (top) 0.3 — 2keV Chandra image with annular sectors from 
region P (labeled counterclockwise from annular sector 1) and 1.5 GHz 
VLA radio contours for NGC 4782(3C278) overlaid. The radio map has 
6.0 X 6.0arcsec resolution with contour levels of 1.5, 7 (cyan), 15, 25, 30, 
35, 40 (green) mJy/beam, showing the radio jets and lobes (green) embedded 
in a faint, large-scale radio halo (cyan), (bottom) Mean surface brightness in 
the 0.3 — 1 keV band as a function of sector number for region P with mean 
radius r = 27kpc. Errors are lcr. The deficit in X-ray emission in annular 
sector 4, coincident with the eastern radio lobe, indicates the presence of an 
X-ray cavity. For comparison, the mean surface brightness measurements for 
annular sectors 2 — 6 at mean radius r = 48 kpc (annulus P„ j ) outside the east- 
em lobe, i.e. 2.2 ±0.2, 2.6 ±0.2, 2.2 ±0.2, 2.0 ±0.2, 2.0 ±0.2, respectively, 
in units of 10~ 4 counts cm -2 pixel -1 , agree within their lcr uncertainties and 
thus show no decrement in X-ray emission between region 4 outside the radio 
lobe and neighboring regions. 1 pixel = 1.968" X 1.968". 

~ 12.6 kpc (40") before flaring into the western radio lobe. 
The western lobe appears distorted into an ellipsoidal shape, 
bending gradually to the northwest, with the radio contours 
on its southern side compressed relative to those on its north- 
ern side. The radio jets and lobes are embedded in a faint, 
large-scale radio halo extending towards the northeast. 

4. 1 . The Eastern X-ray Cavity 

We use the individual annular sectors from the upper left 
panel of Figure[3]to search for X-ray cavities associated with 
the radio lobes. In the annulus (P) extending 17 < r < 37 kpc 
from the center of NGC 4782 and shown in the top panel of 
Fig. |U annular sector 1 contains the western radio lobe, while 
the eastern radio lobe almost entirely fills the annular sector 
4. In the bottom panel of Figure @] we plot the mean sur- 
face brightness in the 0.3-1 keV energy band for each annular 
sector of region P as a function of sector number. For com- 
parison, we also list the mean surface brightness for annular 
sectors 2-6 at larger radii (37 < r < 59 kpc) from region P n \ 
outside the eastern radio lobe, i.e. annular sector 4 and the two 



neighboring annular sectors on either side. We find a 29^}|% 
decrease (lcr uncertainties) in the mean surface brightness in 
sector 4 of annulus P, at the location of the eastern radio lobe, 
while there is no analogous decrease in mean surface bright- 
ness in annular sector 4 of region P„i outside the eastern lobe, 
compared to its neighboring annular sectors. This indicates 
an X-ray cavity coincident with the eastern radio lobe. We do 
not observe a statistically significant decrease in the surface 
brightness (~ 1 %, lcr uncertainties) in the annular sector 
1 of annulus P, that contains the western radio lobe. 

The expected reduction in X-ray surface brightness is small 
for cavities whose inclination angle is large, so that X-ray cav- 
ities are difficult to observe unless the cavity lies close to the 
plane of the sky. We define the surface brightness decrement 
as the ratio of the observed surface brightness in the cavity 
to that expected if there is no cavity, and find the observed 
decrement in the eastern lobe is 0.71^g}| (lcr uncertainties). 
Assuming that the radio plasma in the eastern lobe has emp- 
tied the cavity of X-ray emitting gas, we can use the reduction 
in X-ray surface brightness to constrain the inclination angle 
of the eastern radio lobe with respect to the sky. Modeling the 
expected surface brightness decrement in the 0.3-1 keV band 
by integrating through lines of sight passing through a spheri- 
cal cavity coincident with the eastern radio lobe, as a function 
of the inclination angle of the lobe, we find that the observed 
decrement in the eastern lobe suggests that the eastern lobe is 
inclined at an angle of < 30° with respect to the plane of the 
sky. 

4.2. Jet Emission in X-rays 

In Figure 5 we show the correspondence between the X- 
ray emission and the radio structure of 3C278 in more detail. 
We find X-ray emission associated with each of the two radio 
jets. The top panel shows (as in the top panel of Fig. |4| the 
large scale radio structure from 1 .5 GHz VLA measurements 
superposed on the 0.3-2keV Chandra image for reference. 
In the middle panel, we see that the eastern radio jet bends 
sharply through the 'fan' of X-ray emission extending north- 
east of NGC 4782, with ~ 60% of the 0.3-2keV 'fan' emis- 
sion from the region lying, in projection, within the eastern 
jet. The brightest emission region in the fan (labeled 'blobl' 
in the middle panel of Fig. 5 ) is located near the base of 
the eastern radio plume, and contains ~ 14% (34 ±6 source 
counts within a 3" circular region) of the total 0.3-2 keV fan 
emission. We also find 11 ±3 source counts in the 0.3-2 keV 
band in a 1.5" circular region (blob2) near the bend in the 
eastern jet. If these 'blob' regions are associated with the ra- 
dio jet, we may be observing directly its deceleration and/or 
jet-plume transition. 

In the lower panel of Figure 5 we show a 0.3 - 2 keV Chan- 
dra image of the inner ~ 5 kpc region surrounding the nu- 
cleus of NGC 4782 with radio contours from a 1" x 1" res- 
olution 4.9 GHz VLA map of the same region superposed. 
Two sources, labeled E and W2, are coincident with the ra- 
dio bright ends of both the inner eastern and western jets. We 
estimate the background rates for these two sources using a 
local annulus surrounding each source. We find that these 
sources are soft, with net source counts in the 0.3-2keV 
band of 16.5 ±4.0 and 10.7 ±3.3 (lcr uncertainties) and in the 
2-8keV band of 0.8 ± 1.5 and 0.1 ± 1.0 (lcr uncertainties), 
respectively. Assuming Galactic absorption, these colors are 
consistent with steep photon indices, i.e. T > 2.0. Using the 
observed count rates and this power law spectral model, we 
estimate 1 keV flux densities of < 0.5 nJy and < 0.4 nJy for 
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FIG. 5. — (top) 0.3-2keV Chandra image with 1.5 GHz VLA radio con- 
tours for NGC 4782(3C278) as in Fig. [4] overlaid showing the large scale 
radio structure of 3C178 for reference. The X-ray image has been back- 
ground subtracted, exposure corrected and smoothed with a 1 " Gaussian ker- 
nel, (middle) Zoom-in of the same image as above. 'Blobl(2)' identifies ex- 
tended emission in the X-ray 'fan' possibly associated with the eastern radio 
jet. (bottom) 0.3 — 2 keV Chandra image of the inner ~ 5 kpc of NGC 4782 
with contours from a 1.0 X 1.0 arcsec resolution 4.9 GHz VLA radio map of 
the inner jets of 3C278 overlaid. Radio contour levels are 0.2, 0.4, 0.5, 0.7, 
1.0, 2.0, and 3.0 mjy/beam. E and Wl, W2 label X-ray point sources (knots) 
coincident with the end of the inner eastern jet and the double peaked struc- 
ture of the western radio jet, respectively. The ULX north of NGC 4782's 
nuclear region is also labeled. 



the E and W2 sources. From the VLA radio map, we measure 
peak 4.9 GHz radio fluxes of 4 and 6 mJy for the E and W2 
sources, respectively, and estimate 2-point spectral indices of 
a > 0.9. Although the uncertainties are large due to our low 
statistics, these values are consistent with that expected for 
synchrotron emission from X-ray knots due to the decelera- 



tion of the inner radio jets (see, e.g. Hardcastle et al. 2005b). 
We also find X-ray emission (12 ±3.5 counts when 3 ±0.4 
background counts are expected) extending along the axis of 
the inner western jet, with a knot of X-ray emission (labeled 
Wl in Fig. 5) coincident with a second peak in the 4.9 GHz 
radio emission in this region. 

5. CONSTRAINING THE INTERACTION KINEMATICS 

With the thermodynamic properties of gas in the 
NGC 4782/4783 dumbbell and in the surrounding group de- 
termined, we use these properties to constrain the kinemat- 
ics of the interactions. As shown in Figure Q] the northern 
galaxy NGC 4783 shows a sharp surface brightness discon- 
tinuity ("cold front"; see Vikhlinin et al. 2001; Markevitch 
et al. 2000) on its eastern side. In {35.1l we use the gas proper- 
ties on either side of this surface brightness edge to constrain 
the three-dimensional velocity of NGC 4783. The evolution 
of radio jets and lobes can be strongly influenced by ram pres- 
sure from the surrounding gas, caused by either bulk motions 
in the ISM/IGM or by the motion of their host galaxy thro ugh 
the intragroup medium (see, e.g. Eilek et al. 1984). In {35.21 
and jj5.3l we determine constraints on the ram pressure veloc- 
ity of gas impacting the radio jets and lobes, inferred from 
the thermodynamic properties of the surrounding gas, the jet 
and lobe morphology and the properties of the X-ray cavity 
associated with the eastern radio lobe. 

5.1. Cold Front Analysis Constraints on NGC 4783 's 
Velocity 

Following the method of Vikhlinin et al. (2001), we define 
a bounding ellipse centered on the NGC 4783 nucleus with 
semi-minor (semi-major) axes of 4". 5 (8". 8) and position an- 
gle 304°, that traces the surface brightness edge with the mi- 
nor axis oriented along NGC 4783's direction of motion in 
the plane of the sky. We then construct the surface brightness 
profile, shown in Figure|6l from elliptical annuli concentric to 
the bounding ellipse and constrained to lie in an angular sec- 
tor between 94° and 186° (measured counterclockwise from 
west). This sector was chosen to exclude the 'horn' of emis- 
sion extending northwest from NGC 4783 and emission from 
the bridge south of NGC 4783. In Figure|6l the radius r is the 
weighted mean for each annular region measured from the 
center of NGC 4783 outward along NGC 4783's projected di- 
rection of motion. 

We fit the surface brightness profile across the edge by in- 
tegrating simple spherically symmetric density models along 
the line of sight as a function of the projected distance r. 
We assume a power law model for the density distribution, 
« g ai = n iD (r/ r edge ) ai , for galaxy gas inside the edge (r < r edge ). 
Outside the edge, we assume that the group gas density is well 
represented b y the po wer law fit to the IGM surface brightness 
profile from {33.3.21 i.e. «igm = ^Tigm/ 15.45 kpc) ^ with 
j3 = 0.35 from our previous analysis and hcm measured from 
the center of NGC 4782. The discontinuity (J 2 ) in the X-ray 
emission at the position of the surface brightness edge (r e d ge ), 
is then 

j2 = A i"»in (3) 

An 2 { J 

where n m (n out ) are the gas densities and A; n (A out ) are the 
X-ray emissivities for the galaxy (IGM) gas inside (outside) 
the surface brightness edge, respectively. We then allow the 
edge position (r e d ge ), gas density inside the edge (n; n ), power 
law index inside the edge («i), and the IGM normalization 
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Mach number 

FIG. 6. — (top) 0.3 — 2keV surface brightness profile from elliptical annuli 
concentric to the edge in NGC 4783 and constrained to lie in the angular sec- 
tor between 94° and 1 86° (measured counterclockwise from west) as a func- 
tion of the distance from the nucleus of NGC 4783. The vertical dashed line 
denotes the best fit position of the edge (r ct i sc = 1.76kpc). The solid line de- 
notes the density model fit to the surface brightness for the cavity-uncorrected 
data (Model A) assuming «„ a i oc (r/redge) -Q!1 for galaxy gas inside the edge 

('" < '"edge) and «igm = ^(nGM/15.45kpc) _3/3 , with j3 = 0.35 and hgm the 
distance from the center of NGC 4782, for r > r c( ] oc . (bottom) Pressure ra- 
tio between the undisturbed IGM and pressure inside the edge versus Mach 
number of the galaxy from Vikhlinin et al. (2001). The solid horizontal line 
denotes the pressure ratio derived from the model fit shown in the left panel. 
The dot-dashed horizontal lines denote the uncertainty in the pressure ratio, 
due to the 90% CL uncertainties in the galaxy and IGM temperatures and 
the IGM abundance, and due to variation of the galaxy metallicity over the 
0.3- 1.2Zq range. 

(K) to vary. We find r e( j ge = 1.76kpc, at\ = -1.05, and the 
discontinuity across the edge J = 12.4. This fit (Model A) is 
shown as the solid line in the top panel of Figure [6] 

The velocity of the galaxy is inferred from the pressure ra- 
tio po/pi between the undisturbed (free streaming) IGM and 
that just inside the cold front in the galaxy, that we use as a 
proxy for the IGM pressure at the stagnation point. We use 
the temperature of the IGM gas (kT = 1 .54+.Q ^ keV, abun- 
dance 0.46^q [4Z0), determined from the spectral model fit 



to angular sector 4 in region P, (the region lying directly east 
of the leading edge of NGC 4783), as representative of the 
free streaming IGM in front of NGC 4783. We find a den- 
sity discontinuity across the edge of n m /n oM = 10 ±3 and a 
pressure ratio of po/pi = 3.2±j g, where the uncertainties re- 
flect both the uncertainties in temperature and abundance in 
our spectral fits. We infer that NGC 4783 is moving at Mach 
I.4+07 relative to the IGM, or v = 870^° km s" 1 relative to 
the group IGM. The upper limits on the pressure ratio (5.9) 
and velocity (Mach 1.9, v ~ 1140kms) occur for the case of 
low abundance (0.3Zq), high temperature (0.61 keV) galaxy 
gas in the galaxy's outer gas halo, coupled with low tem- 
perature (~ 1.32keV) IGM gas. The lower limit on the ve- 
locity (Mach~ 0.7, v ~ 470 km s" 1 ) requires high metallic- 
ity (1.2Z ), low temperature (0.37 keV) galaxy gas. Physi- 
cally, this latter combination of parameters is unlikely to oc- 
cur, since such high metallicities are most often present in the 
central cores of elliptical galaxies rather than at the contact 
discontinuity between the outer gas halo of the galaxy and the 
IGM. 

We note that the above hydrodynamic determination of the 
total velocity of NGC 4783 is independent of any possible 
motion of the other interacting elliptical galaxy (NGC 4782) 
relative to the group gas. If, however, we assume that the 
measured relative radial velocity between NGC 4783 and 
NGC 4782 (-623 km s -1 ) is representative of the relative ra- 
dial velocity of NGC 4783 with respect to the group gas (i.e. 
NGC 4782 is the dominant group elliptical galaxy and nearly 
at rest with respect to the group IGM), then we can con- 
strain the three dimensional motion of NGC 4783 through 
the group. This assumption is supported by the dynamical 
measurements of group galaxies by Quintana et al. (1996) 
that place NGC 4782 at the systemic radial velocity of the 
group and is consistent with the X-ray surface brightness pro- 
file of the IGM gas, that is well fit by a spherically symmet- 
ric /3-model centered on NGC 4782 (See 33.3.21 and Jones 
et al. 2007). For a relative radial velocity v r = -623 km s -1 , 
our analysis implies a transverse velocity for NGC 4783 of 
600 (< 960) km s" 1 and an inclination angle C = 46° (> 33°) 
with respect to the plane of the sky. In this model, the motion 
of NGC 4783 must be at least transonic (Mach ^ 0.9 for an 
IGM gas temperature of 1 .54+" keV), since the total veloc- 
ity of the galaxy must equal or exceed its radial velocity with 
respect to the IGM. 

This purely hydrodynamical determination of the veloc- 
ity and inclination angle with respect to the group IGM 
(and NGC 4782) is in good agreement with the binary or- 
bit parameters derived by Madejsky & Bien (1993), who 
find a relative velocity between NGC 4783 and NGC 4782 
of ~ 1000 km s" 1 and inclination angle for the plane of the 
orbit with respect to the plane of the sky of ~ 40° from 
numerical simulations of the tidal distortions of the stellar 
isophotes and velocities expected from the close encounter of 
two equal ~5x 10 u Mq galaxies. The binary orbital parame- 
ters from Borne et al. (1988) favored a lower relative velocity 
(~ 700km s" 1 ) and thus larger inclination angles (£ ~ 65°). 
While our analysis is also formally consistent with this set 
of orbital parameters, the sharpness of the observed surface 
brightness edge for NGC 4783, makes such a large inclina- 
tion angle unlikely (Mazzotta et al. 2001). 

Our analysis for the region east of NGC 478 3 's leading 
edge is complicated by the presence of the radio lobe and 
associated X-ray cavity. We assumed that the apparent co- 
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incidence of the leading edge of NGC 4783 with the outer 
contours of the radio lobe (see Figure @]l is a projection ef- 
fect, such that the pressure difference between the galaxy gas 
inside the cold front and that in the free-streaming IGM is 
due solely to ram pressure caused by NGC 4783's motion 
through the IGM. We make a simple estimate of the maxi- 
mal effect of the cavity on the surface brightness profile in 
front of NGC 4783's edge by repeating the analysis using the 
X-ray decrement to rescale the observed surface brightness 
for radii 2 < r < 20kpc (which overlap the X-ray cavity). 
This estimate for the maximal effect of the cavity on our anal- 
ysis (Model B) reduces the inferred velocity of NGC 4783 
and its upper uncertainty limit (in the steeply rising tran- 
sonic/supersonic regime of Fig. |6| by < 10%, while the lower 
uncertainty limit is lowered by ~ 24%. 

We also investigate whether NGC 4783 will ultimately 
merge with NGC 4782 by using the two body approximation 
to compute the critical mass M a required to lie within the 
radius of the orbit (r ~ 40". 6) for NGC 4783 to be bound 
to NGC 4782. We find a critical mass of (1.1 ±0.8) x 
1O 12 M for Model A, and lower bound of l.O^g x 1O 12 M 
for Model B, respectively. Assuming that the group centered 
on NGC 4782 is in hydrostatic equilibrium and approximately 
isothermal, we can estimate the mass in the central region of 
the group using the /3-model fit to the gas density for the group 
gas. Using the mean IGM temperature of 1 .4keV, we estimate 
a total mass of 7 x lO n M within the central 12.7 kpc (40".6) 
of the group. If a gas temperature of 2 keV is more represen- 
tative of this central region, as suggested by the lower right 
panel of Fig. [3] the total mass is 1O 12 M . After correct- 
ing for differences in the assumed cosmologies, these mass 
estimates are in good agreement with total mass estimates 
for NGC 4782 ((7 ±2) x 10 n M Q ) obtained from the simu- 
lations of Borne et al. (1988), and also with the dynamical 
mass estimates (6 x 1O U M within r ~ 3 kpc) of Madejsky 
et al. (1993). Within the uncertainties, these masses are likely 
sufficient for NGC 4783 to be bound to NGC 4782. 

5.2. Kinematical Constraints from Bending the Eastern 
Radio Jet 

The morphology and evolution of radio jets offer another 
sensitive probe of galaxy and IGM motions. Assuming that 
ram pressure from the surrounding gas dominates, Euler's 
equation relates the ram pressure p ext Vg, where p ext is the gas 
density of the surrounding gas and v g can be interpreted as 
either the relative velocity of the host galaxy through the am- 
bient IGM or the velocity of local bulk flows impacting the 
jets, to the velocity (yf) and density {pj) of the radio plasma 
in the jet, the width of the jet (hj) and the radius of curvature 
R of the bend, i.e. 

pA = p^ (4) 

R hj 

(see, e.g., Eilek et al. 1984; Hardcastle et al. 2005a ). 
The measured minimum (equipartition) pressures of the ra- 
dio plasma in the eastern and western jets are 1.1 x 10" 11 and 
1 .2 x 10" 11 ergs cm" 3 , respectively, and are a factor ~ 3 higher 
than the minimum pressure (4.4 x 10~ 12 ergs cm" 3 ) measured 
in each of the radio lobes. We use the minimum pressure in 
the eastern jet to estimate the minimum density pj = 3/? m j n /c 2 
of plasma in the jet (c is the speed of light), and determine 
the ram pressure velocity. Models for jet deceleration in FRI 
sources find that the flow within the jet is always transonic, 
with relativistic Mach numbers inside the jet of 1 . 1 - 2, corre- 



sponding to velocities of (0.3-0.7)c for light, relativistic jets 
(Bicknell 1994; Laing & Bridle 2002). We adopt the mean 
velocity (0.5 ± 0.2)c as a reasonable estimate of the jet ve- 
locity in 3C278 (NGC 4782), and interpret the observed lack 
of significant beaming between the eastern and western jets 
as further evidence that the jets (and lobes) lie close to the 
plane of the sky. At the bend in the eastern jet, the jet width 
is ~ 6".4 and the radius of curvature of the bend is ~ 16". 
If the eastern jet is passing through the group IGM behind 
the 'fan', we use the /3-model fit to the group IGM around 
NGC 4782 (from €DT2l and Jones et al. 2007) to determine 
the external gas electron density at the location of the bend to 
be 7^Q5 6 x 10~ 3 cm" 3 . From equation @] we find that a veloc- 
ity of v g ~ 1 60 ± 70 km s" 1 would be sufficient to produce the 
observed eastern jet bend. 

The second possibility is that the eastern jet passes through 
the disrupted ISM in the 'fan region' (F). Using the lower 
bound on the density of the galaxy gas in the fan of 9^{J 2 x 
10" 3 cm" 3 (see Table[3]l, derived from our spectral models and 
the assumption of uniform filling, for the external gas density 
acting on the jet in equation |4] we find an upper limit on the 
velocity of v g < 140^™kms _l for ram pressure bending of 
the jet. These velocities are modest, ~ 20-60% of the sound 
speed c s in the cool (0.56 keV) galaxy gas in the fan. Thus, 
while they could be partly due to the motion of NGC 4782 
dragging the jet through the surrounding gas, they could also 
be due to bulk motions in the gas. Madejsky (1992) found ev- 
idence of increased stellar velocity dispersion in NGC 4782 
near the position of the jet bend, evidence that stellar mass in 
that region has been highly disturbed by the tidal encounter 
with NGC 4783. Tidal and pressure forces from that en- 
counter could have also induced bulk motions in the ISM gas 
at the location of the bend with comparable or greater veloci- 
ties, i.e. ~ 140 km s" 1 , capable of bending the eastern jet. 

Since we do not have a direct measurement of the jet ve- 
locity and density of plasma inside the jet, it is useful to 
consider the behavior of our solutions as these parameters 
vary. In typical FRI sources initially light, highly relativis- 
tic (0.8 -0.9c) jets are thought to be decelerated to transonic 
velocities (0.3 -0.7c) by the entrainment of external matter, 
with density contrast Pj/p ex t growing from a few xlO" 6 to 
10" 4 by the end of the deceleration region (Bicknell 1994; 
Laing et al. 1999; Laing & Bridle 2002). For the more pow- 
erful FRI sources, like 3C278 where the onset of flaring is 
> 2 kpc from the nucleus, the deceleration may extend well 
beyond 10 kpc, such that the density contrast at 5-6 kpc from 
the nucleus (at the bend in the eastern jet) is likely interme- 
diate between these extremes. Assuming a density contrast at 
5 kpc of ~ 2 x 10~ 5 in equation|H taken from the jet decelera- 
tion model for 3C31 (see Fig. 7 in Laing & Bridle 2002), the 
ram pressure velocity needed to bend the eastern jet increases 
to 420 ± 170 km s" 1 , still subsonic or transonic with respect 
to the external medium and so not unreasonable in a strongly 
interacting galaxy pair. 

If the jet is much heavier and still fast, jet bending would 
require NGC 4782 to be moving supersonically with respect 
to the group IGM, which is not seen, or that highly supersonic 
fluid velocities impact the jets from the southeast as well as in 
the west, which is also unlikely. If the jets are light, but very 
non-relativistic (e.g. ~ 0.03c as in Colina & Borne 1993), 
then the velocities required to bend the jets would be very low 
(~ 10 km s" 1 ). Such jets would bend too easily and one would 
expect a narrow-angle tail (NAT) morphology rather than the 
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observed mildly bent structure. Finally if the jet plasma is 
heavy, but slow, it would also have to be hot (kT > 20 keV) 
such that the region containing the eastern lobe would still 
appear as an X-ray surface brightness decrement. However, if 
that were the case, the jet would be far from equipartition, and 
the lobe highly overpressured with respect to the IGM. Thus it 
is most likely that the jet plasma is light and relativistic, and is 
bent by subsonic bulk gas flows and/or motion of NGC 4782. 

5.3. Constraints from the Radio Lobes 

Properties of X-ray cavities (bubbles) associated with ra- 
dio lobes place important constraints on the galaxy and/or gas 
velocities in this system. For a cavity filled with relativistic 
plasma, the enthalpy of the cavity is 4pV where p is the pres- 
sure of the hot gas evacuated from the cavity and V is the cav- 
ity's volume. Using the IGM density distribution from 33.3.21 
and assuming the cavity lies in the plane of the sky, we find 
an IGM electron density of 1 .6 x 10~* cm" 3 at r = 73"(23 kpc), 
the radial distance from the nucleus of NGC 4782 to the 
center of the eastern lobe. For an IGM gas temperature of 
1.54keV (from the lower right panel of Fig. |3), the pres- 
sure of the hot gas evacuated from the expanding bubble was 
~7.5 x 10~ 12 ergs cm" 3 . Note that this is less than a factor of 
two larger than the 4.4 x 10" 12 ergs cm" 3 measured minimum 
pressure inside the radio lobe. For a spherical volume of ra- 
dius 34"(10.7kpc), the enthalpy carried in the eastern cavity 
is ~ 4.4 x 10 57 ergs. If the cavity is inclined at 30° with re- 
spect to the plane of the sky (the maximum angle of inclina- 
tion suggested by the X-ray decrement), the density, pressure 
and bubble enthalpy each decrease by < 10%. 

We estimate the age of the cavity in two ways, using the 
geometry of the observed eastern cavity (see, e.g. Birzan 
et al. 2004; Dunn et al. 2005). First we obtain a lower 
bound on the age by assuming the bubble expands at the sound 
speed of the ambient medium (~ 640 km s" 1 ), which is likely 
most appropriate for young bubbles. Second, we use the time 
for the bubble to rise buoyantly, which is likely appropriate 
when the cavity is more evolved. We find a cavity age of 
35-54Myr, within a factor 1.4-3.7 of the time since percen- 
ter passage of NGC 4783 past NGC 4782 (Borne et al. 1988; 
Madejsky & Bien 1993). Thus the episode of AGN activity re- 
sponsible for these features may have been influenced by the 
close encounter of the two galaxies. Using these timescales 
and the enthalpy (4pV) of the eastern cavity, we find a kinetic 
luminosity of ~ 2.6-4 x 10 42 ergss _1 . This is more than 
two orders of magnitude greater than the upper limit on the 
0.5- lOkeV X-ray luminosity for the AGN, estimated from 
the spectral fits in 33.21 but still < 4% of the Bondi accretion 
power (assuming a 10% efficiency) theoretically available to 
the AGN. 

Assuming no intrinsic sidedness, such that the enthalpy of 
the eastern and western lobes are equal, we can use the non- 
spherical geometry of the western lobe to constrain the exter- 
nal pressure exerted on the bubble. The projected geometry of 
the western lobe is an ellipse with semi-major (semi-minor) 
axes of 32" (18"), respectively. Taking the three-dimensional 
geometry of the western cavity to be either an ellipsoid of ro- 
tation about the semi-major axis (prolate) or about the semi- 
minor axis (oblate), we can use the predicted enthalpy of 
the bubble (4.4 x 10 57 ergs s" 1 ) to infer the external pressure 
against which the bubble had to expand. The pressures in- 
ferred from the bubble enthalpy are 2.7 x 10" 11 ergs cm" 3 and 
1 .5 x 10" 1 1 ergs cm" 3 for the prolate and oblate geometries, re- 
spectively. Using the IGM density distribution and mean tem- 



perature (~ 2.5 keV) for gas to the west of NGC 4782 from 
33.3.21 we measure a mean thermal pressure for IGM gas at 
the location of the western cavity of 1.4 x 10~ n ergs cm" 3 . 
The difference between the thermal pressure of the IGM and 
that predicted by the enthalpy carried in the bubble is inter- 
preted as the ram pressure needed to distort the bubble ge- 
ometry. For the prolate geometry this implies a ram pres- 
sure velocity of ~ 600km s" 1 (~ 0.7c s for 2.5 keV gas), a fac- 
tor ~ 3 higher (depending on the assumed jet density), than 
the ram pressure velocity inferred from the eastern jet. How- 
ever, the oblate (pancake) geometry may be more likely, since 
one might expect ram pressure from the southwest to cause 
the bubble to be compressed along the southern rim. For an 
oblate geometry, the predicted ram pressure velocity is lower, 
~ 170 km s" 1 (Q.2c s for 2.5 keV gas). Although the ram pres- 
sure velocity, determined from the lobe analysis, is likely un- 
certain by factors of a few, due to the uncertainties in the lobe 
geometry and IGM gas temperatures, it is intriguing (and grat- 
ifying) that the ram pressure velocity for NGC 4782 required 
to bend the eastern jet (~ 100- 200 km s" 1 ) and western lobe 
(~ 170 km s" 1 assuming oblate geometry) are similar, and so 
likely result from the same dynamical processes. 

6. DISCUSSION 

The interaction picture of NGC 4782 and NGC 4783 that 
emerges from our analysis is the following: The IGM gas tem- 
perature and X-ray surface brightness distribution show that 
NGC 4782 and NGC 4783 reside in a large, massive galaxy 
group (LGG 316) with elliptical galaxy NGC 4782 likely the 
dominant group galaxy near the center of LGG 316's deep, 
dark matter potential. The lack of a sharp X-ray surface 
brightness edge to the south of NGC 4782 supports the in- 
terpretation that NGC 4782 is the dominant group galaxy 
and is not moving rapidly with respect to the center of 
the group. In contrast, analysis of the sharp X-ray surface 
brightness edge on the eastern (leading) side of the ellipti- 
cal galaxy NGC 4783, coupled with its narrow, trailing ram- 
pressure stripped tail, shows NGC 4783 infalling supersoni- 
cally through the group IGM. As NGC 4783 passes superson- 
ically from southwest to east past NGC 4782, tidal interac- 
tions distort the spatial and velocity distributions of stars and 
gas in and between the galaxies. In addition, the compres- 
sion wave produced in the IGM by NGC 4783's high speed 
passage through the group gas impacts NGC 4782's ISM, 
displacing it to the east relative to the stellar isophotes and 
creating the 'fan'. These bulk gas velocities, perhaps aided 
by a small tidally-induced oscillation of NGC 4782 relative 
to the center of the group's gravitational potential, provide 
the ram-pressure needed to bend the radio jets. Our analysis 
demonstrates that understanding the hydrodynamical interac- 
tions between dumbbell galaxies and their surrounding gas 
environment, as revealed by X-ray observations, plays a crit- 
ical role in constraining the interaction kinematics of these 
systems and in modeling their evolution in groups, as well as 
in richer cluster environments. 

7. CONCLUSION 

In this paper we have used the results of a 49.3 ks Chan- 
dra X-ray observation of the dumbbell galaxies NGC 4782 
and NGC 4783 to measure the thermodynamic properties of 
hot gas in and surrounding the interacting system and have 
compared our results with new 1.5 and 4.9 GHz VLA maps 
of the dumbbell. We then used these properties to constrain 
the kinematics of the interaction and models for the bending 
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of the radio jets associated with host galaxy NGC 4782. We 
find the following: 

• In the interaction region close to the two galaxies, the 
X-ray emission is highly asymmetric, showing features 
characteristic of both tidal and gas-dynamical stripping. 
NGC 4783 possesses a sharp leading surface-brightness 
edge 1.76kpc east of its center and a ~ 15kpc tail ex- 
tending to the west, characteristic of ram pressure strip- 
ping, and fixing the direction of motion of NGC 4783 
in the plane of the sky to be from west to east. 

• The temperature, density, and pressure ratios between 
gas inside NGC 4783's leading edge and undisturbed 
IGM gas in front of NGC 4783 indicate that NGC 4783 
is moving at velocity ~ 810^ km s" 1 (Mach 1 .4^7) 
with respect to the group gas. Assuming a radial ve- 
locity of -623 km s _1 ( NGC 4782 radially at rest with 
respect to the IGM), we find an inclination angle for 
NGC 4783's motion of 46°(> 33°) with respect to the 
plane of the sky. This is in agreement with the binary 
orbital parameters of Madejsky & Bien (1993) found 
by fitting tidal distortions in the stellar isophotes and 
velocities. 

• A tidal bridge of X-ray emission extends along the 
line between the centers of NGC 4782 and NGC 4783, 
and a fan of X-ray emission extends to the east of 
NGC 4782. Hot gas in the central regions of both galax- 
ies is well fit by single temperature APEC models with 
kT ~ 0.4-0.5keV, while the tail and bridge tempera- 
tures are higher (kT ~ 0.68 keV). The 'fan' to the east 
of NGC 4782 is composed of cool (0.56 keV) galaxy 
gas. The brightest X-ray emission in the fan is near the 
base of the eastern radio lobe. 

• Outside the interaction region, the X-ray surface bright- 
ness of the group IGM is approximately spherically 
symmetric and nearly isothermal with mean temper- 
ature — 1.4±0.4keV. The IGM surface brightness 
is well fit by a power law density distribution n e = 
2.5+31 x 10 _3 (r/15.45kpc)- 3 ' 3 cm" 3 with = 0.35, cen- 
tered on NGC 4782. However, the data suggest that 
close to the interaction region, at about the radii of the 
radio lobes (17 < r < 37kpc), the mean gas temperature 
to the west of the dumbbell is higher (2.5 I^q keV) 
than the gas temperature to the east (1.75+q jgkeV). 
We find the group X-ray luminosity and total mass 
within a radius of 100 kpc of NGC 4782's center to be 
logLx = 41.99^) ergs s" 1 and (5.5 ± 1.5) x 1O 12 M , 
respectively. Using the M-T relation from Finoguenov 
et al. (2001), we estimate a total mass (at r^oo) for the 
group of 4.8^2 2 x 10 13 M Q , consistent with the large 
(447 km s -1 ) galaxy velocity dispersion measured for 
the group as a whole. X-ray morphology and gas 
temperatures at r > 37 kpc favor the interaction pic- 
ture in which the elliptical galaxy NGC 4783 is in- 
falling for the first time into a single massive group with 
NGC 4782 nearly at rest at the center of the group po- 
tential. 



• There is a 29+J3 % ( 1 cr uncertainties) reduction in the X- 
ray surface brightness in the region coincident with the 
eastern radio lobe of NGC 4782, indicating the pres- 
ence of an X-ray cavity. The observed X-ray surface 
brightness decrement suggests that the inclination an- 
gle of the eastern lobe is < 30° with respect to the plane 
of the sky. Assuming relativistic plasma fills the lobe, 
the enthalpy carried in the cavity is <~ 4.4 x 10 57 ergs. 
The age of the eastern cavity is <~ 35 - 54 Myr, and its 
kinetic power is ~ 2.6-4 x 10 42 ergs s" 1 . 

• The eastern radio jet of 3C278 bends sharply to the 
northeast at ~ 5 - 6 kpc from the nucleus of NGC 4782, 
inflating shortly thereafter into a roughly spherical lobe 
of radius ~ 10.7 kpc, coincident with the X-ray cav- 
ity. The western jet extends straight ~ 12.6 kpc be- 
fore flaring into an ellipsoidal lobe that bends grad- 
ually to the northwest. The minimum radio pres- 
sures in the eastern and western jets are 1.1 x 10 -11 
and 1.2 x 10~ n ergscm" 3 , respectively, and 4.4 x 
10~ 12 ergs cm" 3 in each of the lobes. These equiparti- 
tion pressures are within factors ~ 2-3 of the pressures 
of the surrounding hot IGM gas. 

• X-ray knots are found at the end of the eastern inner 
radio jet and in the western inner radio jet, coincident 
with peaks in the 4.9 GHz radio emission, and are likely 
associated with deceleration of the radio plasma. 

• A ULX with L x ~ 3 x 10 39 ergs s" 1 is found near 
NGC 4782's nucleus. 

• Assuming a light, mildly relativistic jet, ram pressure 
velocities of ~ 100 - 200 km s" 1 are sufficient to pro- 
duce the observed bending in the eastern jet. Assuming 
no intrinsic jet sidedness and an oblate 3-dimensional 
geometry for the western lobe, we find the ram pres- 
sure velocity of the IGM relative to the western lobe to 
be ~ 170 km s -1 with a factor of order ~ 2 uncertainty. 
These velocities may reflect bulk motions established 
in the gas during the galaxies' high velocity encounter, 
with a possible additional small contribution caused by 
the oscillation of NGC 4782 with respect to the group 
potential, also induced by the tidal interaction of the 
galaxy pair. 
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TABLE 1 

VLA OBSERVATIONS USED IN THIS PAPER 



Program ID 


VLA configuration 


Frequency (GHz) 


Date 


Time on source (min) 


AC 104 


C 




1.46, 1.51 


1984 May 03 


11 


AC 104 


c 




4.84, 4.89 


1984 May 03 


147 


AC 104 


A 




4.84, 4.89 


i no A r\,... o/; 
iyo4 Dec ZO 


211 


AC 104 


B 




1.46, 1.51 


1985 May 13 


43 


AC 104 


B 




4.84, 4.89 


1985 May 13 


221 


AH343 


D 




4.84, 4.89 


1989 Nov 17 


9 


Note. — 


All observations were made at two observing frequencies, quoted in this table. 


The effective frequency, the mean 


of these two frequencies, is used 


in the text. 


















TABLE 2 










Spectral Analysis Regions for NGC 4782 and NGC 4783 




Region 


Shape 




Center 


Dimensions 


Orientation 








R.A. Dec. 












(J2000.0) 


(arcsec) 


(deg) 


H 

Hr 


circular 




12 54 36.6-12 33 28.2 


7 






circular 




12 54 37.6 -12 34 36.8 


7 






circular 




12 54 33.9-12 34 36.8 


7 






circular 




12 54 36.3-12 34 46.2 


7 






circular 




12 54 35.1 -12 34 46.2 


7 




T 


rectangular 




12 54 34.8-12 33 36.3 


36,17 


333 


T bg 


rectangular 




12 54 36.6-12 34 37.3 


36,17 


333 


B 


rectangular 




12 54 36.4 -12 33 49.3 


20, 14 


113 


B bg 


rectangular 




12 54 35.1 -12 34 24.5 


20, 14 


113 


F 


rectangular 




12 54 37.3-12 33 59.4 


25, 18 


294 


Fbg 


rectangular 




12 54 36.5 -12 34 29.2 


25, 18 


294 


N 


circular 




12 54 35.7-12 34 06.6 


2.5 




Nbg 


annular 




12 54 35.7-12 34 06.6 


3, 11 




A 


annular 




12 54 35.7-12 34 06.6 


3, 11 




Abg 


annular 




12 54 35.7-12 34 06.6 


12, 25 


195-340 


Note. — 


- Dimensions specified are radii for circular res 


dons, (inner, outer) radii for annular regions, and (length, width) for rectangular regions. Local 


background regions are specified by the subscript bg. Orientation angles are measured counterclockwise from west for the rectangular region's major (length) 


axis and the annular sector A(, g (see Fig. [2}. 
















TABLE 3 








Spectral Fits for Elliptical Galaxies NGC 4782 and NGC 4783 




Region 




Source 




m 


X 2 /dof 






counts 




keV 




H 




195 ± 15 






5.8/8 


T 




290 ±26 






23.6/19 


B 




208 ± 19 




0.68&» 


13.3/10 


F 




146 ±20 




0.56!; 


2.7/10 


A 




213 ±20 






10.9/13 



NOTE. — Col. 2 lists the net source counts in the 0.3 — 8 keV range for regions H, T, A and 0.3-2 keV for regions B and F. The spectral models are absorbed 
APEC models with an additional fixed 5 keV bremsstrahlung component to model the contribution from unresolved LMXBs. The hydrogen absorbing column 
is fixed at the Galactic value («h = 3.58 X 10 20 cm~ 2 ; Dickey & Lockman 1990) for all model components in each region and metal abundances for the APEC 
thermal plasma components are fixed at A = 0.5 Zq . The temperatures are unchanged when the abundance is varied between 0.3 — 1 -2Zq . Errors correspond to 
90% confidence limits. 
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TABLE 4 

X-ray Luminosities for Regions in NGC 4782/NGC 4783 



Region 


£gas 


L gas 


^LMXB 


^LMXB 




(0.5-2keV) 


(2-10keV) 


(0.5-2keV) 


(2-10keV) 


H 


5.4 


0.1 


2.5 


3.8 


T 


9.4 


0.3 


1.3 


2.0 


B 


7.3 


0.3 


1.4 


2.1 


F 


5.0 


0.1 


0.9 


1.4 


A 


3.4 


0.03 


4.7 


6.9 



NOTE. — Intrinsic X-ray luminosities in units of 10 39 ergss -1 for the galaxy emission regions shown in Fig. f2]and spectral models listed in Tableland 
assumed luminosity distance of 66.7 Mpc. The quoted values for the 0.5-2 (2—10) keV luminosities, that assume gas metallicity O.5Z0, change by < 4%(17%) 
when the abundance is varied between 0.3 — 1 .2Zq . 



TABLE 5 

Electron Densities, Thermal Pressures and Gas Masses for Regions in NGC 4782/NGC 4783 



Region 


n e 


P 


M g 


tc 


(1) 


(10 -2 cm- 3 ) 
(2) 


(10~ u ergs cm -3 ) 
(3) 


(10 7 M Q ) 
(4) 


(Gyr) 
(5) 


H 

T 
B 
F 

A 


9 9+0.5 

9+"' 2 
0.9±0.3 


3 3 +LS 

9 9+0.6 

3 3 +M 

l. 3 f 


9 q+06 
+1'7 

4 yAi 
^■ z -1.4 

5 3 +L2 

*-*-1.5 


0.34$! 
0.86«| 

-0.22 

0.89*™ 

77+0-54 
u -"-0.37 



NOTE. — Col.(l): region identifier, col. (2): electron density, col. (3): thermal gas pressure, col. (4): gas mass, col. (5): gas cooling time. Spectral models are 
given in Tablef3]for the spectral regions defined in Tableff] Values are listed for uniform filling and zero inclination angle with respect to the plane of the sky. 
Uncertainties reflect both the 90% CL uncertainty in the gas temperature and the uncertainty from varying the gas metal abundance between 0.3 — 1.2Zq in the 
spectral models. 



TABLE 6 
Nuclear Region Spectral Fits 



Region 


Source 
counts 


m 

keV 


r 


X 2 /dof 


N 
N 
N 


193 ±14.5 
193 ±14.5 
193 ±14.5 


0.58 

43+ - 17 


3.44 
1-7/ 


26/8 
21/8 
7.0/7 


NOTE. — Col. 2 lists the net source counts in the 0.3 — 5keV energy band for region N. The spectral models are absorbed APEC, powerlaw, or APEC + 
powerlaw with an additional 5 keV bremsstrahlung component to model the contribution from unresolved LMXBs. The hydrogen absorbing column is fixed at 
the Galactic value («h = 3.58 X 10 20 cm~ 2 ; Dickey & Lockman 1990) for all model components in each region and metal abundances for the APEC thermal 
plasma components are fixed at A = 0.5 Zq . Errors correspond to 90% confidence limits. / denotes a fixed parameter. 






TABLE 7 
Group Properties 






Row 


Property 




Value 


Units 


(1) 
(2) 
(3) 
(4) 
(5) 
(6) 
(7) 
(8) 
(9) 
(10) 


Cv 
^500 

Mscaletooo 
P 
"IGM 

A 

logLx(lOOkpc) 
logLxteoo) 




447 ±70 

0.66 
7 x 10 13 
(4.8 + | |) x 10 13 

0.35 
(2.5^- 3 ) x 10- 3 
1.4 ±0.4 
0.3-0.5 

41 99+° 13 

42 73+°- 13 
^ z -' J -0.05 


km s~' 
Mpc 

M Q 
M Q 

cirr 3 
keV 

Z Q 

ergss 1 
ergs s~' 



NOTE. — Row (1) LOS velocity dispersion, row (2) rsoo from Eq. 1, Osmond & Ponman (2004), row (3) dynamical mass, row(4) mass from M — T 
scaling relation for groups (Finoquenov et al. (2001)), row (5) and (6) IGM density model parameters (n e = ;iigm('"/ 1 5.45 kpc) -3 ' 3 ), row (7) and (8) IGM mean 
temperature and metal abundance, row (9) measured total X-ray luminosity within r < lOOkpc, row (10) extrapolated total X-ray luminosity within r^oQ. 
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TABLE 8 
IGM Region Geometry 



Annulus 


Radii 


Angular Extent 


Annular Sectors 




^in> ^out 


vstart* ^stop 




(1) 


(2) 


(3) 


(4) 


P s 


15, 54 


185, 25 


6-9,1 


P 


54,118 


345, 345 


1-9 


P„i 


118, 187 


345, 225 


1-6 


P„2 


187, 236 


65, 145 


3.4 


P„3 


236, 305 


65, 145 


3,4 



NOTE. — Col. (1) annulus region identifier, col. (2) inner radius, outer radius of the annulus in arcsec, col. (3) beginning angle, ending angle for the annulus 
in degrees measured counterclockwise from west, col. (4) 40° annular sectors numbered counter clockwise from 345° (see Fig. [5] The annuli are concentric 
with centers at 12 ,1 54 m 35.8\-12 34'7.3" (J2000.0) 



TABLE 9 

Gas Temperatures and Abundances in the IGM 



Region 


Source 
counts 


kT 
keV 


A 

z e 


^/dof 


P s 
P 

Pnl 
P„2 
Pn3 


911 ±36 
2870 ±73 
2269 ± 67 
416±38 
766 ±52 


1 62 +0 - 31 
*■ -0.22 

1.87±0.28 

1 21 +0M 
-Q-2Q 
1 37+0.38 
1 - J '-0.22 


0.55™ 
0.463,4 
0.28:$?? 
o ivw-ii 

-018 

1 5^ iS 
u - 1J -0.09 


47.5/49 
176.7/171 
121/133 
40.3/58 
115/100 



NOTE. — Col. (2) lists the net source counts in the 0.3-8 keV energy band. All regions used an absorbed APEC spectral model for the IGM with hydrogen 
absorption fixed at the Galactic value, «h = 3.58 X 10 20 cnT 3 . 



